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A modified version of the Intelligent Driver Model was used to simulate traffic in the district 
of Afogados, in the city of Recife, Brazil, with the objective to verify whether the complexity 
of the underlying street grid, with multiple lane streets, crossings, and semaphores, is capable of 
exhibiting the effect of critical density: appearance of a maximum in the vehicle flux versus density 
curve. Numerical simulations demonstrate that this effect indeed is observed on individual avenues, 
while the phase offset among the avenues results in damping of this effect for the region as a whole. 
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I. INTRODUCTION 

Over the past decades various models have been pro- 
posed for simulating traffic flow (see e.g. [l|-[6| and ref- 
erences therein), which range on the phenomenon treat- 
ment detail level (microscopic, mesoscopic and macro- 
scopic models), are based on individual events or global 
temporal behavior (discrete and/or stochastic), imple- 
menting real time or off-line simulations. 

Among these, the so called "car- following" models play 
an important role, as they capture much of the realism 
of the phenomenon, with a high level of detail. In par- 
ticular, the "Intelligent Driver Model" (IDM) proposed 
by Treiber and collaborators 0, [1| has attracted much 
attention, as it is capable of reproducing diverse realis- 
tic effects. This model is based on a system of coupled 
non-linear differential equations, where motion of each 
vehicle depends on the position and speed of the adja- 
cent vehicles. The solution of this system of equations is 
only possible through numerical methods. 

However, none of these models has up to date shown 
the ability to reproduce the empirically observed effect 
of a maximum in the vehicle flux versus density curve, 
which represents a fundamental phenomenon from the 
point of view of urban traffic planning. 

In this work we implement a modified version of the 
IDM model, by assuming a truncated Gaussian distri- 
bution for the "desired" individual vehicle speed values, 
on top of the real street map of the suburb of Afogados 
in the city of Recife, Brazil, where we also implement 
real traffic rules. It is found that the composite effect 
of intermittent traffic lights, multiple merging lanes, and 
random "desired routes" of individual vehicles does pro- 
duce a maximum (albeit weak) in the flux versus den- 
sity curve when a single avenue is considered as part of 
a larger composite system, whereas this effect vanes off 
when the measurements are made on the district system 
as a whole. 
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To this end, we have developed a simulator with a 
"Graphical User Interface" (GUI), using animation to 
facilitate the correction of eventual problems, visualize 
the operation and the simultaneous interactions, and aid 
in a presentation which gives credibility to the model. 
Most importantly, this approach enables one to interac- 
tively adapt the model variables, which greatly aids the 
realization of numerical experiments. 

In the following section we give a brief overview of 
the necessary aspects of the traffic theory and the IDM 
model, in the subsequent section we describe the numer- 
ical experiment and present the results, and finally, we 
draw the conclusions. 



II. MODELING URBAN TRAFFIC 
A. Basic concepts 

There are many different aspects of vehicle properties, 
driver behavior, and the surrounding traffic infrastruc- 
ture that may be taken into account when modeling ur- 
ban traffic. For example, the vehicle weight has direct 
impact on the ability of changing the vehicle speed |{|. 
The movement of a vehicle a in a given lane of an avenue 
may be described by considering the vehicle size /q; , ltS 
position in relation to the lane beginning x a , speed v a 
and the acceleration a a . In addition, individual driver 
characteristics influence the unique form in the vehicle 
conduction [9|, [l0[ , where the driver perception and the 
reaction time are the key elements for the efficient vehicle 
operation [9]. Furthermore, individual driver factors in 
the physiologic context may be taken into account. 

Generally speaking, urban traffic represents a highly 
complex non-linear phenomenon, directly related with 
the vehicle and driver properties, the interactions be- 
tween the vehicles, the properties and configuration of 
the traffic infrastructure (streets, lanes, traffic lights etc.) 
and the quantity of vehicles that enter and exit the ob- 
served region per unit time [l2j . 

The traffic flow q is defined p], [l2|, [l3| as the number of 
vehicles N which pass through a detector during a given 
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FIG. 1: Consecutive vehicles in a lane. 
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where in the case of multiple lanes, the composite vehicle 
flow on L lanes is simply given by 
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Depending on the vehicle density and the street capacity, 
traffic flow may be classified as belonging to one of the 
three basic regimes: free flow, synchronized traffic flow 
(when the street capacity is approached) , and "wide mov- 
ing jam" (when the street capacity is exhausted) pi [lij. 

Vehicles in a given lane are ordered due to the impos- 
sibility of overtaking within the lane |l| , the first in will 
be the first out, unless lane switching is permitted. The 
distance S a between the vehicle a and the vehicle a — 1 
immediately ahead, is given by the distance from the 
back bumper of the vehicle in front to the front bumper 
of the following vehicle, represented by 
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where x a represents the vehicle a relative position within 
the avenue, and l a is the vehicle length, as presented 
schematically in Fig. [TJ Taking the instantaneous speed 
of all the vehicles in the region at a given instant, the 
average speed v s is calculated [l| as the arithmetic mean 
of these values 
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Finally, the vehicle density p is defined |l| as the fraction 
of the space occupied by all the vehicles 
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where K is the lane length. 

Empirical studies of traffic flow indicate a determin- 
istic relationship between density p. flow q and the av- 
erage vehicle speed v s p], 0, [l2|, Q~3|, [l5|. This relation 
is fundamental to the traffic flow theory, expressed by 
the equation q = pv s and graphically represented by the 
Fig. [2j where knowledge of two of these three variables 
permits determining the third. However, the functional 
relationship q(p) remains unknown. 




FIG. 2: Fundamental diagram relating the density p to the 
flow q. The flow capacity q cap is reached at the critical density 
p c . The average vehicle speed v s is found as the inclination 
of the straight line that passes through the origin. 

It is seen from Fig. [2] that in the free traffic region, at 
low density values, the flow increases in a linear fashion 
with the increase of density. The congested traffic begins 
when critical density is reached, which defines the flow 
capacity. From this point on, the flow begins to decrease 
with the density increase, corresponding to the regime of 
congested traffic. 

At present there is no consensus in the literature as to 
the functional form of traffic flow versus density curve. 
Several approaches have been proposed up to date which 
end up extracting different functional forms from the 
same data sets. Also, there seems to be no general con- 
sent up to date as to the relevance of individual com- 
ponents of various existing models for the description of 
diverse effects of interest for the understanding of the 
urban vehicular traffic phenomenon. 



B. Intelligent Driver Model 

The "Intelligent Driver Model" (IDM) is governed 
by a system of non-linear differential equations for the 
individual vehicles indexed by a, given by 



dVq 
dt 



= a 



(6) 



where a is maximum (or intrinsic) acceleration, S is the 
parameter controlling the rate at which the desired speed 
is approached, vo is the "desired" vehicle speed, S a is the 
distance to the preceding vehicle a — 1 



S*(v a ,Av a ) = 5o 



v a Av a 
max Iv a H , 



2Vab 



(7) 



is the "safety distance" between the vehicles a and a — 1, 
5o is the minimum security distance, T is the breaking 
reaction time [16], Av a is the difference between the ve- 
locities of the two vehicles, and b is the deceleration pa- 
rameter governing the influence of the velocity difference 
on the safety distance. 
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Beside an abundance of parameters (perhaps, over- 
abundance for the current purpose), this model contains 
intuitively attractive components for the description of 
vehicular behavior, and has been shown to display rather 
realistic behavior in diverse situations (see [7|, |8| and 
citing references). The solution of the model does not 
present physically unacceptable situations (such as neg- 
ative velocities, vehicle overlapping etc.), and stable so- 
lutions are attainable already with the Euler's method, 
while application of the Runge-Kutta method yields ac- 
curate results for describing effects such as e.g. traffic 
density waves propagation. 

On the other hand, this model yields realistic effects 
only in the presence of obstacles or conflicting situations 
such as confluent lanes, whereas in the case of free lane 
traffic it demonstrates trivial behavior (all vehicles tran- 
sit in parallel with the same speed) , without limit for the 
street capacity. It turns out that using a truncated Gaus- 
sian distribution [13, [H| for maximum desired speeds of 
individual vehicles, is sufficient to create a traffic con- 
gestion in an avenue without obstacles [19], where street 
capacity becomes apparent through saturation of the flow 
curve. However, the maximum of the flow curve at the 
critical density was not observed with this generalization 
of the original IDM model [l9j . 

In order to investigate whether this phenomenon may 
be brought about by the additional confusion stemming 
from the complexity of the environment (multiple lanes, 
vehicle types, traffic rules etc.), in what follows we im- 
plement a numerical simulation implementing such a sit- 
uation. 



III. NUMERICAL EXPERIMENT 

The area of study is located in the Recife Metropolitan 
Region (RMR), where a network of urban highways con- 
nects municipal districts and suburbs. According to Per- 
nambuco State Traffic Department (DETRAN) data, in 
the year of 1990 the RMR was traversed by some 250000 
vehicles, which has grown to over 700000 in 2007, where 
over 400000 of these are found in the city of Recife alone. 
This violent vehicular density increase has been causing 
heavy traffic jams in several parts of the city over the 
past years, augmenting the attention of the authorities 
and researchers. The study region itself represents one of 
the two principal transport corridors between the south 
metropolitan region and Recife downtown. 

As already mentioned, the IDM model by itself does 
not yield results demonstrating critical behavior on the 
flux versus density scatter plot, and the generalization of 
IDM which considers a distribution of desired velocities 
leads only to saturation of the flux versus density curve, 
with no sign of diminishing flux when the density is fur- 
ther increased. In order to verify whether the complexity 
of an urban traffic system is able to reproduce such an 
effect, in the current numerical experiment we have im- 
plemented various components encountered in the real 



life situation. In particular, we consider the following 
basic elements: 

• The avenues are considered in terms of length, 
number of lanes, crossings, and traffic lights; 

• Vehicles of different size and maximum acceleration 
are considered (trucks, busses and cars), with in- 
dividual maximum and minimum "desired" speed. 
Upon entry, each vehicle is (randomly) assigned a 
"desired" itinerary. 

• Traffic rules are implemented in each avenue and 
lane, maneuvering is implemented in accord with 
the existing signalization, considering the possibil- 
ity of changing lanes and overtaking slower vehicles, 
while avoiding collisions. 




FIG. 3: Street map of the Afogados suburb with indicated 
traffic flow directions. 

We present in Fig. [3] the street map of the Afogados sub- 
urb with six entries and seven exits, where the directions 
of the traffic flow are indicated. All the avenue crossings 
were visited and examined physically, in order to take 
into account all the possible maneuvering options of a 
vehicle in a given lane of a given avenue, which attempts 
to follow a given route, while honoring the existing traffic 
rules. 

The traffic lights were modeled as elements that belong 
to the avenues, such that the vehicles detect the traffic 
light presence, and behave in accord with the stop or go 
state (in fact, this was implemented by inserting a virtual 
still vehicle at the crossroad when the traffic light is red, 
such that approaching vehicles sense its presence, and 
break according to the IDM rules). 

As already mentioned, the vehicles entering the region 
are attributed (random) desired routes. In the attempt 
to follow this route, a vehicle chooses permitted maneu- 
vers, attempting to choose a lane that optimizes access 
to a determined crossroad, and realizes lanes changing 
and overtaking in the attempt to reach the destination. 
To this end we implement a origin — destination route 



4 



table taking into account all the possible combinations of 
the entry and exit avenues that comply with the traffic 
rules. While attempting to follow the desired route, the 
vehicle may also change lanes due to advantageous con- 
ditions (e.g. when directly in front of it there is a slower 
vehicle, and the adjacent lane is empty), while honoring 
traffic rules. This basic setup creates a rather realistic 
behavior of individual vehicles, and it seems to capture 
the collective behavior of the vehicle fleet, which is ob- 
served on the graphical interface, as shown in Fig. 2] . 




FIG. 4: A snapshot of the graphic interface of the simulator 
during the simulation. 

In order to achieve a real time simulation on the cur- 
rent scale that permits realistic traffic behavior, we have 
opted for the Euler algorithm with a 0.05s time incre- 
ment, and a full scene refresh at the rate of two frames per 
second. While a higher order Runge-Kutta method yields 
higher precision results, the computational demand on 
this simulation scale takes it out of the real time ob- 
servation range on the current hardware, and we have 
verified that the overall behavior of the traffic system is 
not affected by the precision gain. 

For the model parameters we adopt the values a = 
1.5ra/s 2 , b = 2.0ra/s 2 , T = 1.2s, S = 2.0m, and 5 = 4, 
while the values of the desired maximum speed for each 
vehicle are drawn from a normal distribution truncated 
at Vmi n and V ma x^ depending on the vehicle type, and in 
accordance with the local traffic rules, with parameters 
given in Tab. HI 



TABLE I: Parameters of the normal distribution truncated 
at Vmin and V ma x, used to draw the desired maximum speed 
for different vehicle types (the values are given in km/h). 



Vehicle type Vmax 


Vmin 




a 


cars 80 


40 


60 


20 


buses 50 


30 


40 


5 


trucks 40 


30 


35 


2.5 



Occupation (%) 



Occupation (%) 



The simulation is performed by inserting vehicles into 
all the incoming avenues (and lanes) at the edge of the 



FIG. 5: Simulation results for the avenue 00. 



simulation zone. The insertion rate is gradually increased 
from 0.125 cars per second per lane (8s delay between in- 
coming vehicles), to continuous influx (no delay between 
inserting the vehicles), with successive decrements of 0.1s 
for the delay between incoming vehicles. At each rate (al- 
together eighty values), ample equilibration time of 1800s 
(30mm) is used to attain equilibrium, after which the 
vehicle count and individual speeds are recorded in each 
avenue segment bounded by street crossings and conflu- 
ence points. This process was repeated 100 times in order 
to establish the average values. 



IV. RESULTS 

Results of the simulation corresponding to the avenue 
labeled "00" in Fig. [3] are presented in Fig. [5j where it is 
observed that the flux versus occupation density scatter 
plot does display signs of a maximum (and thus existence 
of critical density), while the average speed versus occu- 
pation demonstrates large fluctuations in the low density 
regime, which diminish with density increase. 

In order to alleviate the effect of large fluctuations and 
pinpoint the critical occupation density, in Fig. [6] we dis- 
play the same data as in Fig. [5j together with the results 
of averaging the data in 2% bins, and polynomial regres- 
sion which indicates a critical occupation of 41.63%, with 
a maximum flow of 0.91998 cars/ s (or 3312 cars/h). 

Results integrated over the entire simulation region are 
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FIG. 6: Polynomial regression of the third order for simulation 
data, avenue 00. 



displayed in Fig. [7J which does not reveal a maximum in 
the flux versus occupation density scatter plot. Also, 
variation of the average velocity is pronounced over the 
whole observed range of density occupation values. This 
findings may be attributed to the fact (observed on the 
graphical interface during the simulation) that the entire 
region never becomes entirely congested, rather, different 
regions take turns in demonstrating free flow, and traffic 
jam states. 



Occupation (%) 



suburb in Recife, Brazil, in order to verify whether the 
additional complexity brought about by the synchronized 
collective behavior of vehicles that attempt to follow in- 
dividual routes, is capable of producing the critical occu- 
pation density effect with a maximum of the flux versus 
density curve. This behavior is empirically observed, and 
represents a fundamental effect from the point of view of 
urban traffic planning, but has not been reproduced up 
to date by existing mathematical and numerical models. 

The modification of the IDM model implemented in 
this work consists of several ingredients. In particular, 
we consider three types of vehicles with different charac- 
teristics (length, maximum acceleration, and individual 
desired speed), multiple lanes with traffic rules (merg- 
ing traffic, semaphores etc.), and predefined individual 
routes for every vehicle. 

It turns out that different parts of the observed re- 
gion go intermittently through phases of free flow and 
congested traffic behavior. If only a single avenue is ob- 
served, the composite effect of interaction with the neigh- 
boring regions is reflected in behavior that may be iden- 
tified as critical density (albeit weak in the current sim- 
ulation). However, the phase offset between the regions 
(while a given region is congested, a neighboring region 
displays free flow, and vice versa) leads to canceling out 
of individual congestion effects, and the maximum in the 
flux versus density curve is not observed. Also, the av- 
erage velocity as a function of occupation displays het- 
eroscedastic behavior, with large fluctuations at low den- 
sities that diminish as the density is increased, whereas 
large average velocity fluctuations are observed for the 
entire occupation density range for the suburb as a whole. 

We may conclude that the critical density phenomenon 
should be regarded as a local effect, brought about by the 
interplay of a given region with the neighboring regions, 
as an exit point for one region represents an entry point 
for another. To the best of our knowledge, the current 
work represents the first report in the literature regarding 
a mathematical/numerical model capable of exhibiting 
this effect. 

Further studies should be made as to the contribution 
of each of the implemented model components to the ob- 
served critical density effect, and whether it may be en- 
hanced by inclusion of some others. It should also be 
investigated whether other microscopic models are capa- 
ble of demonstrating this phenomenon in a similar setup. 
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In this work we implement a modified version of the 
Intelligent Driver Model on top of a real street map of a 
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